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Abstract 

We investigate the finite-temperature behavior of the Yukawa model in which Nf 
fermions are coupled with a scalar field <f> in the limit Nf — > oo. Close to the chiral 
transition the model shows a crossover between mean-field behavior (observed for Nf = 
oo) and Ising behavior (observed for any finite Nf). We show that this crossover is 
universal and related to that observed in the weakly-coupled cf) 4 theory. It corresponds 
to the renormalization-group flow from the unstable Gaussian fixed point to the stable 
Ising fixed point. This equivalence allows us to use results obtained in field theory and 
in medium-range spin models to compute Yukawa correlation functions in the crossover 
regime. 
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1 Introduction 



The finite-temperature transition in QCD has been extensively studied in the last twenty 
years and is becoming increasingly important because of the recent experimental progress in 
the physics of ultrarelativistic heavy-ion collisions. Some general features of the transition, 
which is associated with the restoration of chiral symmetry, can be studied in dimensionally- 
reduced three-dimensional models [HI21- However, a detailed understanding requires a direct 
analysis in QCD. Being the phenomenon intrinsically nonperturbative, our present knowl- 
edge comes mainly from numerical simulations [31 II]- Due to the many technical difficulties — 
finite-size effects, proper inclusion of fermions, etc. — results are not yet conclusive and thus 
it is worthwhile to study simplified models that show the same basic features but are signif- 
icantly simpler. In this paper we shall consider a Yukawa model in which Nf fermions are 
coupled with a scalar field through a Yukawa interaction. The action of the model in d + 1 
dimensions is 

N f 

S = DN f /d d+1 x(i(^) 2 + |0 2 + A^ + ^ /d^x^ + ^ + M)^, (1) 

where tr 7 2 = D (D = 2 d / 2 if d is even, D = 2^ d+1 ^ 2 if d is odd), the integration is over 
M. d x [0,T _1 ], and A > to ensure the stability of the quartic potential. Along the thermal 
direction we take periodic boundary conditions for the bosonic field (f) and antiperiodic ones 
for the fermionic fields ipj. The theory must be properly regularized. We shall consider a 
sharp-cutoff regularization, restricting the momentum integrations in the spatial directions 
to p < A. However, the discussion presented here can be extended without difficulty to any 
other regularization that mantains at least a remnant of chiral symmetry. 

In the limit Nf — > oo this model can be solved analytically and one finds that there is 
a range of parameters in which it shows a transition analogous to that observed in QCD 
jSUHJ. It separates a low-temperature phase in which chiral symmetry is broken from a high- 
temperature phase in which chiral symmetry is restored. For Nf = oo this transition shows 
mean-field behavior, in contrast with general arguments that predict the transition to belong 
to the Ising universality class. This apparent contradiction was explained in Ref. [7j where, 
by means of scaling arguments, it was shown that the width of the Ising critical region scales 
as a power of 1/Nf, so that only mean-field behavior can be observed in the limit Nf = oo. 
An analogous behavior was observed in a generalized O(N) a model in Ref. for finite 
values of iV the transition was expected to be in the Ising universality class, while the N = oo 
solution predicted mean-field behavior. In Ref. we performed a detailed calculation of the 
1/N corrections, explaining the observed behavior in terms of a critical-region suppression. 
The analytic technique discussed in Ref. [9 j can be applied to model (JH). It allows us to obtain 
an analytic description of the crossover from mean-field to Ising behavior that occurs when 
Nf is large and to extend the discussion of Ref. [7j to the case M^O. More importantly, 
we are able to show that the phenomenon is universal. In field-theoretical terms, it can 
be characterized as a crossover between two fixed points: the Gaussian fixed point and the 
Ising fixed point. This implies that quantitative predictions for model ([TJ can be obtained in 
completely different settings. One can use field theory and compute the crossover functions 
by resumming the perturbative series [TUl HH H2 EE] • Alternatively, one can use the fact 
that the field-theoretical crossover is equivalent to the critical crossover that occurs in models 
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with medium-range interactions f3J El El El EI- This allows one to use the wealth of 
results available for these spin systems jTH El El El El El- In this case the interaction 
range R is essentially equivalent to a power of N, N ~ R d . Finally, we should note that the 
phenomenon is quite general and occurs in any situation in which there is a crossover from the 
Gaussian fixed point to a nonclassical stable fixed point. For instance, similar considerations 
have been recently presented for finite-temperature QCD in some very specific limit [2Ti] . 

The paper is organized as follows. In Sec. 121 we review the behavior in the limit Nf = 
oo. In Sec. 01 we consider the 1/Nf fluctuations and determine the effective theory of the 
excitations that are responsible for the Ising behavior at the critical point. These modes are 
described by an effective weakly-coupled 4> 4 Hamiltonian. In Sec. 14.11 we present a general 
discussion of the critical crossover limit. These considerations are applied to the Yukawa 
model in Sec. 14.21 and in Sec. 14.31 We determine the relevant scaling variables and show 
how to compute the crossover behavior of the correlation functions. Finally, in Sec. we 
present our conclusions. In the appendix we discuss the relations among medium-ranged spin 
models, field theory, and the Yukawa model. A short summary of this work was presented 
in Ref. [2T]- 



2 Behavior for Nf = oo 

The solution of the model for Nf = oo is quite standard. We briefly summarize here the 
main steps, following the presentation of Ref. [H] . As a first step we integrate the fermionic 
fields obtaining an effective action <S e ff[</>] given by 



-DN f S eS [ 



/=i 



(2) 



where 



S cff [0] = /d d + 1 x(^) 2 + ^ 2 + 



4! 



1 

D 



+ M). (3) 



For Nf — ► oo one can expand around the saddle point 
equation 



A 3 
jum H — m 
6 



m 



+ M)TJ2 [ 



d d+1 x tr lo£ 

- 4>, that is determined by the gap 
1 



/ P <A (27r)> 2 + W 2 + ( m + M) 2 ' 



(4) 



where we define the frequencies u n = (2n + 1)ttT, and 

m = glj) JI = fig- 2 A = A#~ 4 
The action corresponding to a saddle-point solution m is: 



S eS (m,M,T) 



li 2 A 4 T 
—m H — -m 
2 4! 



p<A 



d d p 1-2 +u 2 + / m + M \21 



(2 



7T 



v z + < 



(5) 
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where we have added a mass-independent counterterm to regularize the sum j^j. Such a 
quantity has been chosen so that the action for M = m = vanishes. Summations can be 
done analytically [6 j. The gap equation can then be written as 



p(m) = (m + M) Q(m + M,T), 
where the functions p(m) and Q(x,T) are defined by 

p{m) = Jim + — m 3 , 
6 



(6) 



(7) 



G(x,T) 



d d p 



l P <A (2vr) d 
Analogously we can rewrite Eq. (jSJ) as 



;V V+x 2 /t 



+ i 



5 eff (m,M,T) 



fl 2 A 4 

— m H — r-m 
2 4! 



— T 



p<A 



d d p 

(2n) d 



log 1 cosh -^-^ ^ - I — log cos! i 
i 2T / 2T 



P 



(9) 



Using Eqs. (JHJ) and © we can determine the phase diagram of the model. Given /Z and A, 
for each value of T and M we determine the solutions m of the gap equations. In general, we 
find either one solution or three different solutions mo, m + , and m_ with m_ < mo < m + 
(for some specific values of the parameters two of them may coincide). When the solutions 
are more than one, the physical solution is the one with the lowest action S e s(m, M, T). 

Note that Eqs. (jOJ) and © are invariant under the transformations m — ► — m and M —>■ 
—M. Thus, we can limit our study to the case M > 0. There are four different possibilities: 

(a) If /Z > ^(0, 0) = CqA 4 - 1 with 



Co 



2V /2 (rf 



' .2 



n -1 



(10) 



then, for every M > 0, there is only one solution m + > 0; for M = we have m + = 0. 

(b) If < JZ < CoA d_1 , there exists a critical temperature T c (ji). For T > T c (/Z) and any 
M there is only one solution m + > (for M = we have m + = 0). For T < T c (/Z) and 
M < M there are three solutions mo, m+, and m_ with m_ < mo < m + and m + > 
and m_ < 0. For T < T c (/l) and M > M there is only one solution corresponding to 
m + . The physical solution is always m + so that M has no physical meaning. Moreover, 
for T < T c (/Z) and M = 0, m + > 0. The critical temperature can be computed from 



the following relation: 

7Z=£(0,T C ) = T C J2 



p<A 



d d p 1 

(2n) d p 2 + Lul, 



A d - L C - A 



d-l 



d d p 



where ui^ n = (2n + l)nT c and T c (p) 



p<1 (27r) d p[eP^ + iy 
t c (~fi)A. For fl — > 0, we have T c (/Z) 



oo. 
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(c) If —G\ <~p < 0, with 1 

3/2 , 1/2 3A d rfd + iy 1 (12) 

there is a critical mass M such that there are three solutions for M < M, two of 
them coincide for M = M, while for M > M the only solution is m+. The physical 
solution — the one with the lowest action — is always m + so that M has no physical 
meaning. Note that m + > for M = 0. 

(d) For ~p < —Ci there are three solutions for all values of T and M. The relevant solution 
is always m + > 0. 

In case (a) chiral symmetry is never broken, while in cases (c) and (d) chiral symmetry is 
never restored. Thus, the only case of interest — and the only one we shall consider in the 
following — is case (b), in which there is a chirally-symmetric high-temperature phase and a 
low-temperature phase in which chiral symmetry is broken. 

The nature of the transition is easily determined. We expand 

G(x,T) = YQ mn x 2m (T-T c ) n , (13) 



m,n 



where 



9oo 



9oi 



r< 

E 



d d p 



net 



p<A (27r)<V + u; c y 



• p<A (2rrY (p» + ul n f 
d d p e p/Tc 



T C 2 J P<A (27r) d (e^ + i)2^ 



9io 



.TYf ^ - (14) 



and u Ctn = (2n + l)nT c . Since ^oo = G(0,T c ) = \i [Eq. (|TT|) ]. the gap equation © becomes: 

\rr? = -pM + g 01 (T-T c )(m + M)+g 10 (m + M) 3 + --- (15) 
o 

where we have neglected subleading terms inm + M and T — T c . Defining 

u^-^M u ^^n {T -T c ), (16) 

bgio - A bgio - A 



1 C7i is the solution of the equation p(—x) = ]hxiM— >oo MQ(M, 0) with a; = (2C1/A) 1 / 2 . The value x 
corresponds to the position of a maximum of p(m) for ~p — —C\. 
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and taking the limit u^, u t — > at fixed x = u t /u 2 J 3 we obtain the equation of state 

m = u h ^f{x) (17) 
= f( x f + xf{x) + 1. (18) 

Note that the prefactor of T — T c in w t is always positive to ensure the existence of only one 
solution for T > T c . Such an equation is exactly the mean-field equation of state that relates 
magnetization (p, magnetic field h, and reduced temperature t. Indeed, if we consider the 
mean-field Hamiltonian 

t u 

n = h^ + -^ + - v \ (19) 

the stationarity condition gives 

11 

h + tif + -^ = 0, (20) 
o 

which is solved by <p = Ah 1 ^ 3 f(Bt\h\~ 2 ^ 3 ), where f(x) satisfies Eq. IjlHj) . and A and B are 
constants depending on u. This identification also shows that M plays the role of an external 
field, while m ~ is the magnetization. 



3 Effective theory for the zero mode 

In order to perform the 1/Nf calculation, we expand the field <fi around the saddle-point 
solution, 

0(x) = + ^=0(x), (21) 
gVN 

where N = DNf, and </>(x) in Fourier modes: 

0(x d , x d+1 ) = TJ2 e 2 ™ nTx ^ [ ^ 4(p)e Jp - Xd . (22) 



In the following we will refer to the integers n — or more precisely to 2irnT — as frequencies. 
In this way we obtain the following expansion for the effective action: 



S cS [4>n] = N{S cS [4>}-S cS [4>}} 



P<A 

im (Pi) • • • <MpO^ (0 (Pi> ni; • • • ; pz, nj). (23) 

Note the Kronecker 5 on the frequencies that ensures that vertices are nonvanishing only for 
Y2i n i = 0- m particular — this property will be important below — if only one frequency is 
nonvanishing, we have V®(pi, n; P2, 0; • • • ; pz, 0) = 0. 
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The vertices appearing in expansion (|23|) are easily computed. The fermion contribution 
is obtained by considering the one-loop fermionic graphs in theory (0). If we define the free 
fermion propagator 

~ { J2j=i I3P3 ~ Hd+iVn + m 



A F (p, n;m) 



P 2 +UJn + 171 



(24) 



with u n = (2n + 1)ttT, the fermion contribution is 



d d q 

(2n) d 



tr 



( i % 

Y[ A F (q + ^p,; a + J2 n f m + M ) 

i=l j=l j=l 



+ permutations, 



(25) 



where the permutations make the vertex completely symmetric [there are (/ — 1)! terms]. 
Note that the frequencies u n never vanish and thus the vertices have a regular expansion in 
powers of m + M. Vertices vf^ satisfy an important symmetry relation. First, note that 



7 d+ iA F (p, n; m)^ d+1 = -A F (p, -n - 1; -m). 



(26) 



It follows 

aeZ 



m) 



1 = 1 



E 



tr 



II 7d+iAi?(qi; a + &*; m)j d+1 



i=i 



(27) 



■i)'E 



tr 



aG2 



-a — 6,- — 1; — m) 



4 = 1 



tr 



a62 



I J A/ (q,: a - b { ; -m) 



In the second step we used 7^ +1 = 1, while in the last one we redefined a — > —a + 1. This 
relation implies (we write here explicitly the dependence of the vertices on m and M) 2 



Vj l \pi, m; . . . ; pi, n z ; m + M) = (-l) i l// ; (p 1 , -m; . . . ; Pi, -m; -m - M). (28) 

For every / > 4, the vertex is due only to the fermion loops, so that V® = vj l \ For / < 4 
we must also take into account the contribution of the bosonic part of the action, so that 

V (3) (p,n l ;q,n 2 ;r,n 3 ) = Am + V} (3) (p, m; q, n 2 ; r, n 3 ), 
^ (4) (p,ni;q,n2;r,n 3 ;s,n4) = A + V^ 

Finally, for the inverse propagator we have 



A + V f (4) (p, m; q, n 2 ; r, n 3 ; s, n 4 ) 



(29) 



P(p,n) 



p 2 (2ixTn 2 A 2 Tr(21/ 

^ + - 2 -^ + /I+-m 2 + K; ; -p,-n;p,n . 

9 9 2 1 



(30) 



2 If c? is odd, one can repeat the same argument using 7^+2 = Il<=i 7»' ^ snows that vertices with Z legs 
are multiplied by (—1)' if one changes the sign of m + M at fixed momenta and frequencies. 
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Vertices also satisfy the symmetry relation (|28|). while the inverse propagator P(p, n) 
satisfies P(p, n; m, M) = P(p, -n; -m, — M). Finally, note that 1/ (4) (0, 0; 0, 0; 0, 0; 0, 0) is 
positive at the transition. Indeed, one obtains explicitly (note that A > to ensure the 
stability of the quartic potential) 



V*»(0, 0; 0, 0; 0, 0; 0, 0) = A + 6T e £ / > 0. (31) 



' q<A 

It is easy to verify that P(0, 0) vanishes at the transition. Indeed, for m = M = we have 

d d q 1 



P(0,0) = JZ-TJ2 [ 



l q <A (2tt)V + ^ 
= -g 01 (T-T c ) + O(T-T c ) 2 , (32) 

where we used Eqs. (II 3|) and (fTT|) . Thus, the mode with n = is singular. It is exactly 
this singularity that forbids a standard 1/Nf expansion at T = T c and gives rise to the Ising 
behavior. This type of singular behavior is completely analogous to that observed in Ref. [§]. 
The strategy proposed there consists in integrating all the nonsingular modes <fi n , n^0, and 
study the effective theory for the zero mode O - 

Integrating all fields n with we obtain the effective action 



e -S eff [0 o ] 



\{^ ri e-s^\ (33) 

n^0 



with 



4 ff = Vnhmo) + ^ [ £^Mp)Hp)M-p) (34) 

A Jq<A 

, \- T 1 - 1 f d d Pl d d Pl d 

+ gfl^y«<A^-'w (27r) ^S R 

x\/W(p 1 ,...p / )0 o (p 1 )---0o(p/) 

where PT, P, and V™ have an expansion in powers of 1/N. The computation of these 
quantities is quite simple. The contribution of order 1/N k to V"> is obtained by consid- 
ering all /c-loop diagrams contributing to the /-point connected correlation function of 4>q 
and considering only the nonsingular fields (i.e. propagators with n ^ 0) on the internal 
lines. Frequency conservation implies that all tree-level diagrams with more than one vertex 
vanish. 3 Therefore, H = O^ 1 ), P(p) = P(p, 0) + C^iV" 1 ), and V® = V® =0 + C^iV" 1 ) 

(Ki,'=o * s the ver t ex V® with all frequencies set to zero). For the inverse propagator P and 

3 For a tree-level graph, the usual topological arguments give the relation J2 n ( n ~ 2)JV„ = —2, where N n 
is the number of vertices belonging to the graph such that n legs belong to internal lines. Since n > 1 if 
there is more than one vertex, the previous equality requires iVi > 2. But frequency conservation implies 
that Vi vanishes if all frequences but one vanish. Therefore, each nontrivial tree-level diagram vanishes. 
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for the magnetic field H we shall also need the 1/N corrections. We obtain 

H = § + 0(iV~ 2 ) (35) 
P(p) = P( p ,0) + ^M + O(N-i), (36) 



with 



* - fE/ <A (§^(P,«)P(P,,)-' (37) 



p(1) (0) = |W ^k(p ) n)P(p 1 n)- 1 -%n) J P(p l n)4 (38) 

where VJ( P , n) = V^(p, n; — p, — n; 0, 0; . . .). Note that relation (J25j) implies an analogous 
symmetry relation for V®: 

V®( Pl , ...,pi;m, M) = (-l)'y«( Pl , . . . , Pj; -m, — M), (39) 

where we have written explicitly the dependence on m and M. Analogously P(p) and if 
are respectively symmetric and antisymmetric under m,M —> —m, —M. 

In order to obtain the final effective theory we introduce a new field x(p) such that the 
corresponding zero-momentum three-leg vertex vanishes for any value of the parameters. For 
this purpose we write 

a X (p)=TMp) + VNk5(p), (40) 

where a and k are functions to be determined. If we write ai = V"\0, 0; . . . ; 0, 0), k is 
determined by the equation 

(-l) l k(m,M) 1 , 
^ 1 ; | ' J at +3 (m, M) = 0, (41) 

1=0 

where we have written explicitly the dependence on m and M. Now, symmetry (}3*9"j) implies 
also 

^(-l) l [-k(-m,-M)} 1 . 

E „ -<W™, M) = 0, (42) 

i=o L 

so that k(m, M) = —k(—m, —M). Therefore, k has an expansion of the form 

k= k ^ a M\ (43) 

ab,a+b odd 

where the coefficients k a b have a regular expansion in powers of 1/N . The leading behavior 
close to the transition is easily computed: 

k= — + 0{m a M b ,a + b = 3). (44) 

Q-4 
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In terms of \ the effective action can be written as 

S m = N^m x (0)+ l - j i^ x (p)P(p) x (-p) (45) 

+ E L §F ' " & (2l)d < £ p> ) F< " (Pl1 ' ' ' p,)x(Pl) ' ' ' * (p,) ' 

The quantities H, P, and have an expansion in terms of to, M, and 1/N. Explicitly we 
have: 

h 2 k 3 

H = aT' l [H - fcP(O) + — V 3 (0) V 4 (0) + 0{m a M\ a + b = 5)], (46) 

2 6 

P(p) = ^"^(p) - A;V 3 (p) + y!4(p) + 0(m a M\ a + b = 4)], (47) 

F (2 ' +1) ( Pl , . . • , P2m ) = a^T-^^dh, . . . , p 2m ) - kV^( Pl , p 2l+1 , 0) 

+0{m a M b ,a + b = 3)], (48) 

F (2 °( Pl , . . . , p 2i ) = a^T-^WCpx, . . . , P2i ) + 0(m a M b , a + b = 2). (49) 
Up to now we have not defined the parameter a. We will fix it by requiring 

dP(p 



dp 2 



= 1, (50) 

p=0 



for all values of the parameters. The parameter a is a function of to, M, and 1/N. The 
symmetry properties of k and of the vertices imply that a is invariant under m,M — > 
—to, — M. As a consequence, under to, M — > —to, — M, the quantities if, P, and have 
the same symmetry properties as if, P, and V"®. 

In the following we shall need the expansions of H, P(0), and V^(p, — p, 0) close to the 
critical point. For this purpose we will use the expansions 

P(0, 0) » ^to 2 — (T — T c )(? 01 - 3(M + to) Vo, 
\/ 3 (0, 0) w Ato - 6(M + m)g w , 

V 4 (0,0)« A-6<? 10 , (51) 

[K(p, 0) = \/W(p, 0; -p, 0; 0, 0; 0, 0, . . .)] and the relation 

Vf (3) (p,n;-p,-n;0,0;m) = mV^p, n; -p, -n; 0, 0; 0, 0; to) + 0(to 3 ), (52) 

where we have explicitly written the mass dependence of the vertices. We expand H and 
P(0) is powers of 1/N as 

™=fc„ + ^ + 0(iV- 2 ), (53) 

ZZW =p0 + | + O (iV-). (54) 
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By using expansions (|51|) we obtain 

v 3 (o,Q) n/n , i v 3 (o, Q) 3 

« -JIM + -^M(T - T c ) - ftoA(6g 10 + A) M 3 + 0(m a Mb; a + 6 = 5); 
65-10 - A (65-10 - A) 2 

AMT x ^ f d d p 

"" n) 



Po « P(0,0) 



2(6010 - A) f^ J P <A 

0{m a M b ,a + b = 3) 
V 3 (0,0) 2 



2V 4 (0,0) 



-001 (T - T c ) + %loA _ M 2 + 0(m a M\ a + b = 4), 
65-10 - A 

P (1) (0,0) « e + 0(m a M 6 ,a + 6 = 2), (55) 



Pi 



where e is the value of P^'(0, 0) for M = m = 0. Note that several terms that are allowed by 
the symmetry m, M — > — m, — M are missing in these expansions. In the case of ho we used 
the gap equation to eliminate the term proportional to m 3 . This substitution is responsible 
for the appearance of the term linear in M and cancels the terms proportional to m(T — T c ), 
m 2 M, and mM 2 . In the case of hi and po note that the terms proportional to m, and m 2 , 
mM cancel out. Finally, we compute the three-leg vertex. At leading order in 1/N we obtain 

Jf' 3 »(p,-p,o) = k 3 (p.o)-|||v4(p,o) 

_/U/ [65i + V r f (4) (p,0;-p,0;0,0;0,0)] + O(m a M b ,a + 6 = 3). (56) 



65-10 - A 

Note that the term proportional to m is missing as a consequence of relation (|52j) . 



4 The critical crossover limit 

The manipulations presented in the previous section allowed us to compute the effective 
action for the zero mode x(p). F &r from the critical point P(p) 7^ for all momenta 
and thus one can perform a standard 1/iVf expansion. At the critical point instead this 
expansion fails because P(0) = 0. At the critical point, for N 00 the long-distance 
behavior is controlled by the action 



S, 



efi 



d d X 



-(d X f 



u 
4! 



X 



+ 0(N~ 



(57) 



where 



u = ^v (i \o, 0,0,0). 



(58) 
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Here we have used the fact that vertices with an odd number of fields vanish at the critical 
point and the normalization condition (|50|). Moreover, since the critical mode corresponds to 
p = 0, we have performed an expansion in powers of the momenta, keeping only the leading 
term. Since for N -> oo, F (4) (0, 0, 0, 0) = c^T^V^O, 0, 0, 0), inequality flflj implies u > 
at the critical point. Eq. (J57j) is the action of the critical 4 theory which should be studied 
in the weak-coupling limit u — > 0. In this regime the model shows an interesting scaling 
behavior — we name it critical crossover — that describes the crossover between mean-field 
and Ising behavior. This allows us to obtain quantitative predictions for the critical-region 
suppression observed for Nf — > oo. 



4.1 The general theory 

In this section we wish to review some basic results on the critical crossover limit. An 
extensive discussion can be found in Refs. [TOl EH El • Let us first consider the standard </> 4 
theory in d dimensions with d < 4. 



5 m „tM = / err 



(59) 



We assume the theory to be regularized with the introduction of a momentum cutoff. The 
results are however independent of the chosen regularization and one could equally well use 
a lattice regularization. Then, we define a new bosonic field i[>(s) as 

^(s) = uM/PW-^su-VM). (60) 
Formally, the action can be rewritten as 



Scant = J d d S 



(61) 



where 

H EE i?w -(<i+2)/[2(4-<i)] ) ~ _ m -2/(4-d)_ ( q 2) 

Thus, formally, once the action is expressed in terms of ip, the bare parameters appear only 
in the combinations H and f. Then, consider the zero-momentum connected correlation 
function Xn- We have 



Xr 



e j d d r 2 . . . d d r n (y?(0Mr 2 ) . . . ^(r n )) conn 

= u [2d-n(2+d)]/[2(4-d)] J ^ _ _ _ ^ ^( )V>(S 2 ) • • • V(Sn)) C ° nn 

= M [ M -"( 2+ ^/[ 2 ( 4 -^/ n (iJ,f), (63) 

i.e. jj-l^-^P+^/Pt 4 -^]^ i s a scaling function of H and f. Analogously, one can determine 
the scaling behavior of the correlation length £: 

e= 2dY 2 I dW ^(°M r )> = ^^ 2/(4 - d) / € (^,r). (64) 
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The above-reported discussion is valid only at the formal level since we have not taken 
into account the presence of the cutoff that breaks scale invariance. For d < 4 only a 
mass renormalization (a redefinition of the parameter r) is needed in order to take care 
of divergencies. By a proper treatment [10] one can show that there is a function r c (u) 
such that the correlation function \ n and £ satisfy the scaling relations (J53j) and (JMj) with 
t = tu~ 2 ^ 4 ~ d \ t = r — r c (u), replacing f: 

Xn = u N-nM]/[2(4- d ) ]/n(#) f) ^ = u -m^) M ^ t). (65) 

The function r c (u) takes care of the ultraviolet divergent diagrams. In two dimensions, only 
the tadpole is primitively divergent and we have 

r c (u) = —\nu + Ku. (66) 

In d = 3 divergences appear at one and two loops, so that 

A u 2 

r c (u) = u + —— Inu + Ku 2 . (67) 

In both cases the arbitrary constant K can be chosen so that t = corresponds to the critical 
point. The function r c (u) depends on the chosen regularization (the expressions we report 
above correspond to a sharp-cutoff regularization). On the other hand, the scaling functions 
f(H, t) are regularization- independent (universal) once a specific normalization for the fields, 
the coupling constant, and the scaling variables is chosen. They are the crossover functions 
that relate mean-field and Ising behavior. Consider, for instance, the case H = 0. For t 
fixed and u — > we obtain the standard perturbative expansion; thus, t — > oo corresponds 
to the mean-field limit. On the other hand, for t — > at u fixed, Ising behavior is obtained; 
t = is the nonclassical limit. By varying t between and oo one obtains the full universal 
crossover behavior. 

In Ref. jH] we extended these considerations to the general two-dimensional Hamiltonian 
S M [cp] = Jty(0) + iy^[tf(p)+rMpM-p) (68) 

+£^ir / f^"W {2nfs fe v " ,(Pl ' ' ' ' ' p,Mpi) ' ' ' v(p,) ' 



where K(p) =p 2 + 0{p 4 ), V (3) (0,0,0) = 0, and V (4) (0, 0, 0, 0) = 1. The presence of vertices 
with an odd number of legs requires an additional counterterm for the magnetic field. Indeed, 
we showed that it was possible to find functions r c (u) and H c (u) such that for t = r — r c {u) 
(infrared limit), h = H — H c (u), u — > (weak-coupling limit), at fixed t/u, h/u one has 

Xn = u 1 - n f n (h,t), (69) 

where the scaling function f n (x,y) is the same as that computed in the continuum theory. 
In particular, Xn'u n ~ 1 vanishes in the crossover limit if n is odd. The counterterms are 
regularization-dependent. In the continuum theory with a cutoff we have 

f d' 2 p V 3 (p) .... 
u \ u ^ u f d 2 p V 3 (p) 2 
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with V 3 (p) = V (3) (p, -p, 0) and 
A = 



3 1 , 3 1 
D 2 ~ — + — log 



vr ' 8tt ~~° 8tt 2 J p<A (2tt) 2 



d 2 p 



y( 4 )(p,- P ,o,o) 

K(p) 



P 



The nonperturbative constant D 2 was estimated in Ref. |13| : D 2 = —0.0524(2). 



(72) 



4.2 Scaling behavior 

In this section we wish to use the previous results to compute the crossover behavior of 
model (JIJ in 2+1 dimensions. Since u ~ 1/N the relevant scaling variables are 



X h 
X t 



NTc -(N l ' 2 H - H c ) 



a 
NT C 

a 2 



(P(0)-r c ), 



(73) 



where the factors a/T c and a 2 /T c are introduced for convenience. The critical crossover limit 
is obtained by tuning T, M, and N close to the critical point so that Xh and Xt are kept 
constant. The expansions of H and P(0) are reported in Eq. ()55|). The expansions of iif c 
and r c are easily derived. For H c we have 



f(3) 



d 2 p V w (p,-p,0) 



2v / iVy p <A(27r) 2 P(p) 
where all quantities are computed for M = m = 0. Using Eq. (|56J) we have 



ah c0 M 



+ 0(m a M b ,a + b = 3), 



T C VN 

where h c o is a constant. Using Eq. (fTTj) we obtain for r c the expansion 



a 



T C N 



(r In + n) + 0{m a M\ a + b = 2 



where 



r 



a 2 1/ 4 (0,0) 



57T 



aV 4 (0,0) 



1 , 3a 4 1/ 4 (0,0) „ 3 
— hi — — — Do 

Stt 87rT r A 2 8vr 



d 2 p 



T c ^ 4 (p,0) a 2 V4(0,0) 



(2tt) 2 L P(p,0) 



P 



(74) 



(75) 



(76) 



(77) 



(7f 



Note that the three-leg vertex that appears in Eq. (|7T|) does not contribute to this order, 
since it vanishes for m = M = 0. Thus, Eqs. (|75j) can be written as 



x t 



N 



-JIM + a M(T - T c ) + aiM 3 + 

e 



-(7oi(T-T c ) + a 3 M^ + - + 



a 2 M 

— + " 

• ro In iV — ri 



h c0 MVN 



(79) 
(80) 
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where ciq, ai, 0,2, 03, and e are coefficients that can be read from Eq. (|55p. These expansions 
show that 

JIM = -x h N-V\ (81) 

m m e — r In N — r 1 — x t , . 

T-T c = — i ^. 82 

The critical point is specified by the condition xt = Xh = 0. The symmetry under m, M — > 
—m, —M guarantees that the critical point corresponds to M = 0. On the other hand, 1/N 
fluctuations give rise to a shift of the critical temperature. If T C (N) is the finite- N critical 
temperature, we obtain 

Te{N)xTc+ ^^hILzIl . (83 ) 

Note that, beside the expected 1/N correction there is also a In N/N term that is related to 
the nontrivial renormalization. It follows 

T-T e (N) = --p-. (84) 

Note that g i is negative [see Eq. (fTljl] and thus we have x t > for T > T C (N), as expected. 
Using the gap equation we can also derive the behavior of m in the critical crossover limit. 
We obtain 

where m is a function of In iV that satisfies the equation 

i(A - 6g w )ml + (r In iV + r x - e + x t )m + x h = 0. (86) 


For iV — > 00, mo has an expansion in inverse powers of In N, the leading term being 

m ^-^-±- + O(\n- 2 N). (87) 
To m A/ 

Note that m — > as Xh — > 0. 

These results confirm the scaling predictions of Ref. |7]. For the massless theory with 
M — 0, there are two regimes: for N(T — T C (N)) 1 one observes Ising behavior, while 
for N(T — T C (N)) ^> 1 mean-field behavior occurs. If M / the same considerations 
apply, the relevant variable being MN 3 ^ 2 . It is important to note the role played in the 
derivation by the symmetry m,M —>■ —m, —M, that is present because the regularization 
preserves chiral invariance. Even though vertices with an odd number of legs are present, the 
symmetry makes them irrelevant in the crossover limit. Thus, the additional renormalizations 
computed in Ref. |9] do not play any role here. 

The results reported above can be extended to d dimensions for d < 4. Eq. (|52*jl implies 
that the relevant scaling variables are 

x h ~ MN^I^ A - d ^ , Xt~[T- T C {N)}N 2 /^ . (88) 

In d = 3, on the basis of Eq. (|67|). we also predict for T C (N) an expansion of the form 

m /, T N m o&hiiV-|-c ,. 
T C (N) « T c + - + — , (89) 

where a, b, and c are constants that can be computed as in the two-dimensional case. 
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4.3 Correlation functions 



The results reported in Sec. 14.21 allow us to compute the scaling behavior of the correlation 
functions. For instance, we have 



(<j>{^ d ,x d+1 )) = - - + - 7 =( X (x d )) (90) 

9 gVN 

Using Eq. (J55J) with n = 1 and d = 2, we have (x(x d )) = f\{xh,x t ) in the critical crossover 
limit. The background term can be neglected in the crossover limit since 



9 9 



y 3 (Q,Q) 

m Vk(0,0) 



M 6# 



10 



# A - 6pi 



iV~ 3/2 . (91) 



Thus, we obtain 



a 

(^(xd.acen-i)) » —=f x {x h ,x t ) (92) 
OVA' 



The factor 1 / viV is related to the particular normalization of used in (0) and disappears if 
we redefine ip = y/N<p in order to have a canonical kinetic term for cp. The function fi(xh, x t ) 
is the scaling function for the magnetization in the Ising model. For instance, for Xh = 
and xt < (low-temperature phase), we have fi(0,Xf) ~ (— x-tf 1 and fi(0,x t ) ~ (— x t )^ MF 
respectively for |x t | -C 1 and |x t | ^> 1, where /?/ = 1/8 and = 1/2 are the magnetization 
exponents in the Ising and in the Gaussian model. The universality of the crossover allows us 
to compute the scaling functions in any other model in which there exists a crossover between 
the Gaussian and the Ising fixed point. In particular, we can use the results for systems with 
medium-range interactions f3J (see also the appendix). In Ref. jTH] (LBB) the authors 
report (\m\R) versus tR 2 (see their Fig. 9), where m is the magnetization, t the reduced 
temperature, and R the effective interaction range. These results give us (0(x d , Xd+i)) for 
Xh = 0. One only needs to take into account the different normalizations of the fields, of the 
coupling constant, and of the scaling variable. In the crossover limit iV — > oo, T — > T C (N) 
at fixed N(T - T C (N)) we have 

gVN((f>(x d ,x d+1 )) = K 1>hBB (\m\R) LBB (93) 
(^ 2 )lbb = K LBB N[T - T C (N)) . (94) 

The nonuniversal constants K^ BB and K\^ BB are computed in the Appendix. 
It is customary to define an effective exponent f3 e ff(T) as 

/3 eff (T) = [T-T c (N)}^\n(<p(x d ,x d+1 )), (95) 

for M = and T < T C (N). In the crossover limit T -> T C (N), N -> oo at fixed N[T-T C (N)] : 
the exponent f3 e g(T) interpolates between the Ising value j3i = 1/8 and the mean-field 
f3 mf = 1/2. Again, this effective exponent can be derived from the results of Ref. ^H]- The 
curve reported in Fig. 15 of Ref. [IHj gives j3 e g in the Yukawa model once tR 2 is replaced by 
K LBB {T-T C (N)]N. 



16 



The same considerations apply to the connected zero-momentum n-point function Xn'- 



Xn 




d d x 2 ...d d x„(x(0)x(x 2 )... X (x n )) 



N' 



x n (0(O)0(x 2 )...0(x n )) 



n/2-1 



conn 



conn 



(96) 



9 



,n 



For n = 2 the crossover function for Xh = can be obtained from the results of Ref. |18j . 
since g 2 X2 = -^2,lbb(x-R 2 )lbb- The constant ^2,lbb is given in the Appendix. 

One can also use field theory to compute the crossover curves and thus use the results of 
Ref. For instance, in the high-temperature phase, for M = we have in the crossover 
limit 



where F x {t) is reported in Ref. ^3] and Kft, -^2,ft are nonuniversal constants computed in 
the appendix. 

In the discussion presented above we have focused on the case d = 2, but it is immediate 
to generalize all these considerations to the three-dimensional case. For d = 3 the universal 
crossover curves have been computed in Refs. PUJ El E] (field theory) and in Ref. [T9"] 
(medium-range models). These results apply directly to the Yukawa model. 

5 Conclusions 

In this paper we have considered the Yukawa model in the limit Nf — > oo, focusing on 
the crossover between mean-field and Ising behavior. For this purpose we have determined 
the action of the mode that becomes critical at the transition. In the long-distance limit, 
it becomes equivalent to that of a weakly coupled </> 4 theory. This identification allows us 
to use the results available for this model [TO"! |Hj and, in particular, to identify a universal 
critical crossover occurring for Nf — > oo, M — > 0, and T — T C (N) — > at fixed x t and Xh, 
see Eqs. fISljl . (|MJ) . In field-theoretical terms, this behavior represents the crossover 

induced by the flow from the unstable Gaussian fixed point to the stable Ising fixed point. 
Quantitative results for the Yukawa model can be obtained by using the field-theoretical 
results of Refs. jTUl CHI 112] , or the Monte Carlo results available for medium-ranged models 
|18| ITD*] . The necessary nonuniversal renormalization constants can be computed in pertur- 
bation theory. Results for d = 2 are reported in the appendix. 

We should stress that our results are not specific of the chosen regularization, but can 
be extended to other regularizations as well. In particular, the extension to Kogut-Susskind 
fermions, the model considered in Ref. |7j, is essentially straighforward. The Wilson case is 
more involved. Indeed, the absence of chiral symmetry implies that the symmetry relations 
satisfied here by the effective vertices [see Eq. (J28|) ] are no longer valid. In turn, this may 
imply additional mixings as it happens in the generalized Heisenberg model 0. 

Let us note that all calculations presented here refer to the model in infinite spatial 
volume. However, the crossover behavior can also be observed in the finite-size scaling 
limit. The discussion in Sec. 14.11 can be easily extended to this case too. It is trivial to 



g 2 X2 = K 2 ,ftF x {1) 



i = K FT N[T - T C {N)} 



(97) 
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verify that the correct scaling variable is L = Lu l ^ A ~ d \ i.e. L = LiV -1 ^ 4- ^ in the Yukawa 
model. Again, one can use universality and obtain predictions for the Yukawa model from 
the results obtained in other contexts. In particular, one can use the finite-size scaling results 
of Refs. fUE3EI] that refer to medium-range models at the critical point. 

Finally, we should mention that one could also generalize the model and consider fermion 
fields ip a f transforming according to a representation of a group G and a coupling of the form 
ijj fT a ([) a ipf, where T a are the generators of the algebra of G. The discussion is essentially 
unchanged, though in this case one would obtain the vector 4 theory. Field-theory results 
relevant for this case are given in Refs. ^01 113 E2] • 

A Relations among the Yukawa model, medium-range 
models, and field theory 

In this appendix we relate the weakly coupled 4> 4 theory, medium-range models, and the 
Yukawa model for d = 2. For simplicity, we only consider the case H = 0, corresponding to 
M = in the Yukawa model. The field-theory model has been discussed in Sec. 14. 1[ where 
it was shown that the n-point zero momentum connected correlation function Xft> shows a 
scaling behavior of the form 

« n_1 XFT,n = /ft>(*ft) ^ft = (r - r c (u))/u. (98) 

Next, we consider systems with medium-range interactions. Consider a square lattice, Ising 
spins a x at the sites of the lattice, and the Hamiltonian 

We assume 4 that J(x) = 1 for \x\ < R m , J(x) = for |x| > R m . The behavior of these 
models is very similar to that observed in the Yukawa model, R m playing the role of N. For 
any finite R m , the system belongs to the Ising universality class, while for R m = oo all spins 
are coupled together and one obtains mean-field behavior. In Ref. [H] it was shown that 
this model shows a crossover that interpolates between mean-field and Ising behavior. If one 
defines an effective interaction range R by 

tf= E^-y)V(*-y) 
E x , y J{x-y) 

then for R, R m — > oo, t = (T — T C (R))/T C (R) — > at fixed ^mr = tR 2 one has 

R A 3 ™XMR,n = /mR,ti(^Mr)) (101) 

where Xmr,u is the connected zero-momentum n-point correlation function of the fields a. 
In Ref. [12] it was shown that /mr,„(s) and fFT,n( x ) are closely related. Indeed, we have 

fuR,n{x) = ^i.mr^mr/ft^IAmr^), (102) 



4 One can also consider a much more general class of medium-ranged models, see Ref. 
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where /i^mr and Amr are mo del- dependent constants that reflect the arbitrariness in the 
definitions of the fields, of the range R, and of the scaling variable t. The constants can be 
computed using the results of Ref. [T3], Sec. 4. 2. 5 Explicitly, we have for the two-point and 
four-point zero-momentum connected correlation functions: 



XMR,4-R 8 



^MR 



In 



/ 47Tt MR \ 

V 3 J 



+ 8nD 2 + 3 



0(i 



~ 3 ) 

MRJ' 



0(i 



MR) 



In the field-theory model we have instead 

1 



UX2 



1 

t + 87rt 2 



' 8vrt 

In + MD 2 + 3 

3 



+ o{r 



u 2 xa 



TA +o(t- 5 ) 



Comparing we obtain 



A*1,MR 



Mmr — A 



MR 



1 

2" 



(103) 

(104) 
(105) 

(106) 



In Sec. E]we have shown a similar relation for the Yukawa model. If x t = — go\N[T — T C (N)], 
we find for the zero-momentum correlation functions of the field 6 



Xn — 9 Xr, 



N n/2 ~' fy, n {x t ) 



(107) 



and 



fY,n( X ) = ^Y^rf FT, n(^YX). (108) 

In order to compute these constants we compare the one-loop expansions of the two-point 
function in field theory and in the Yukawa model. In the Yukawa model we find 

2 



d d x( X (0) X (x)) 



NT [NT 



a 2 x t 
1 

~2N 



\a 2 x t 
NT 
a 2 x t 



d 2 p V 

WT 2 



( 4 ), 



p,-p,0,0) 



+ 0(x t 



Using the explicit expression for r c we obtain 

8irx t 



X2 



1 a 2 V A (0,0) 



X t 



6TTX7 



hi 



3a 2 y 4 (0, 0) 



+ 8ttD 2 + 3 



+ 0{x t 



For the four-point function we have instead 

V 4 (0,0) 



XaN- 1 



x 



+ 0{x- 



(109) 



(110) 



(111) 



5 Note that the function f x {t) defined in Ref. ^3] refers to correlations of the fields (j> and not of the original 
fields <p. However, relation (4.12) of Ref. [H] shows that in the critical crossover limit Ylx(f°f x ) ~ Ylxi^ ^) ■ 
The same holds for \i- The expression reported here are obtained from those reported in Ref. T!?| by setting 
a,2 = l,a,4 = -2, N = 1, Co = Co = r, and £mr = t + cq. 
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Comparing we obtain 

A i 1 , y = aV 4 (o,0), ^y = ^^Jy A ^ = ^W ^ 

The constants reported in Sec. 14.31 are easily derived: 

K FT = —Xy9oi, K 2: ft = Hi,yH2,yi (113) 

K Xy901 K ^Y^Iy , , 

-"■LBB = 7 , -"-n.LBB = ~ , l 114 J 

A MR /^ljMR^.MR 

where n = 1,2. Note that g i is negative, so that K FT and -Klbb are positive as expected. 
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